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摘 要 : 本 文 以 吉林 省 松花 江 流 域 为 研究 区 ,以 1995 年 .2005 年 和 2015 年 3 期 Landsat TM/ETM 影像 为 数据 源 ,对 


吉林 省 松花 江 流 域 景观 格局 \ 生 态 脆 弱 度 的 时 空 分 布 和 变化 特征 及 驱动 力 进行 研究 ,结论 如 下 :中 流域 不 同时 
候 变 化 存在 差异 ,景观 优势 类 型 以 耕地 和 林地 为 主 ,城乡 用 地 和 草地 稳定 增长 ,林地 


期 气 


地 \ 水 域 呈 减少 趋势 。@) 流 


域 脆弱 度 以 较 低 脆弱 区 和 中 度 脆 弱 区 为 主 , 低 脆弱 区 和 中 度 脆 弱 区 主要 分 布 在 流域 中 游 和 下 游 较 发 达 地 区 ,高 脆 
弱 区 主要 分 布 在 流域 东南 地 区 和 下 游 西北 荒漠 化 及 林 草 地 破坏 区 域 。@) 流域 脆弱 度 总 体 相 对 稳定 , 呈 略 微 上 升 和 


两 极 分 化 的 趋势 ;脆弱 度 与 人 为 活动 密切 相关 ,与 地 形 有 一 定 相 关 性 ， 


问 题 o 


受气 候 变 化 影响 并 不 显著 ,但 仍 能 反映 部 分 
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20 世纪 80 年 代 以 来 ,景观 格局 受到 越 来 越 多 
的 学 者 关注 ,成 为 全 球 环境 问题 的 关键 内 容 和 当前 
国内 外 研究 的 热点 之 一 ”。 景 观 作 为 多 种 栖息 地 
斑 块 的 有 机 结合 体 , 为 人 类 生存 提供 适宜 的 环境 , 具 
有 很 强 的 综合 性 和 鲜明 的 区 域 性 。 景 观 格局 是 指 不 
同形 状 和 大 小 的 斑 块 在 空间 上 的 各 种 排列 组 合 ,是 
景观 异 质 性 的 具体 体现 ,对 生态 过 程 有 重要 意 
义 ,能 够 较 好 地 反应 区 域 生态 系统 的 脆弱 度 ” 。 
目前 ,国内 外 对 于 景观 格局 研究 多 集中 于 LUCC 和 
景观 格局 变化 “生态 安全 与 生态 风险 等 领 
ROI ,脆弱 度 研究 并 不 多 见 , 孙 才 志 等 对 下 辽河 
平原 , 张 月 等 对 艾 比 湖 , 梁 二 敏 等 对 新 疆 玛 纳 斯 河流 
域 绿洲 进行 了 景观 格局 脆弱 度 分 析 ”-” ,这 些 研究 
多 是 对 研究 区 格 网 化 计算 景观 指数 再 进行 插值 ,名 
略 了 景观 太 度 效应 问题 ,而 且 研 究 多 为 脆弱 度 的 格 
局 分 布 和 空间 变化 ,分 析 其 空间 自 相 关 性 ,很 少 结合 
人 为 和 自然 因素 来 考虑 景观 脆弱 度 分 布 格局 形成 与 
变化 的 原因 。 

景观 格局 在 某 种 程度 上 决定 着 生态 与 环境 的 质 
量 , 与 各 生态 过 程 密切 相关 ,对 系统 稳定 性 和 生物 多 
样 性 有 重大 影响 。 景 观 脆弱 度 体 现 了 系统 受到 
外 界 的 干扰 和 脆弱 程度 ,是 评价 区 域 生态 的 有 效 指 
xU 。 河 流 是 人 类 生存 与 繁衍 的 基础 ,流域 生态 系 
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统 的 完好 是 人 类 发 展 和 经 济 运行 的 重要 保证 ,流域 
生态 的 问题 也 是 当前 关注 的 焦点 " ,其 景观 格局 及 
脆弱 度 的 研究 对 于 流域 环境 保护 及 合理 开发 具有 重 
要 意义 ,松花 江 是 我 国 重点 流域 ,也 是 吉林 省 的 母亲 
河 ,20 a 来 气候 及 环境 发 生 显著 变化 ,并 引起 了 草地 
沙化 和 盐 碱 化 等 问题 的 。 本 文 对 1995 年 .2005 年 
和 2015 年 松花 江 流域 吉林 省 部 分 的 LandsatTM 和 
ETM + 影像 数据 进行 解 译 ,得 到 研究 区 这 3 个 时 期 
的 景观 分 布 图 ,通过 建立 相应 的 景观 格局 指数 构造 
生态 脆弱 度 模 型 ,利用 Fragstats ArcGIS 等 软件 , 通 
过 移动 窗口 法 获取 1995—2015 年 松花 江 流域 (吉林 
省 ) 景 观 脆弱 度 时 空 分 布 和 变化 ,并 结合 人 为 和 自 
然 因 素 , 对 驱动 力 进 行 研究 。 研 究 成 果 对 指导 松花 
江 流域 景观 资源 的 合理 开发 和 利用 ,以 及 生态 恢复 
具有 一 定 的 指导 作用 和 应 用 价值 。 


1 研究 区 概况 


松花 江 流域 位 于 119° 52’ ~ 129° 30'E ,41° 42' ~ 
51? 38'N 之 间 , 流 域 面积 5.57 x 10 km2 。 其 中 吉林 
省 松花 江 流 域 面积 约 9. 49 x 10* km? ,范围 见 图 1。 
地 势 的 总 体 趋势 呈 西 北 低 、 东 南 高 。 流 域 西 北部 
于 平原 丘陵 地 貌 , 海拔 在 100 ~200 m, 其 中 北部 
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图 1 研究 区 区 位 图 


Fig.1 Geographical location of the study area 


于 低 山 丘 陵 区 ,南部 属于 长 白山 山系 ,西部 为 松 嫩 平 
原 , 是 松花 江 流域 的 主要 农业 区 域 ,海拔 50 ~ 200 
m, 流域 具有 明显 的 大 陆 性 季风 气候 ,冬季 严寒 漫 
长 ,夏秋 季 降 雨 集 中 ,春季 干燥 多 风 ,年 内 温差 较 大 。 
流域 多 年 的 平均 气温 在 4 %C 左右 ,平均 降水 量 在 
670 mm 左右 ,存在 明显 的 空间 差异 ,呈现 出 由 南 向 
北 递减 的 趋势 。 流 域 矿产 丰富 ,土地 肥沃 。 虽 是 黑 
龙 江 的 支流 ,但 对 东北 地 区 的 人 民生 活 及 农业 生产 
的 意义 都 超过 了 黑龙 江 和 其 他 河流 ,因此 ,流域 的 生 
态 人 研究 非常 重要 ,选取 吉林 省 辖区 部 分 作为 研究 区 。 


2 数据 来 源 与 研究 方法 


2.1 数据 来 源 与 处 理 

本 研究 通过 地 理 空 间 数据 云 网 站 ,选取 吉林 省 
松花 江 流域 1995 年 .2005 年 .2015 年 时 相 为 6—10 
月 的 Landsat TM 和 ETM + 影像 数据 ,每 年 15 景 影 
像 图 ,植被 覆盖 率 较 高 ,影像 包含 的 信息 丰富 , 云 覆 
盖 度 小 ,影像 清晰 容易 判读 ,通过 ENVI 软件 进行 几 
何 校正 ,图像 增强 等 处 理 ,结合 野外 调查 和 相关 资料 
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进行 室内 解 译 ,得 到 3 期 吉林 省 松花 江 流 域 景观 类 
型 分 布 图 。 通 过 矢量 化 获得 流域 道路 和 城市 点 及 水 
系数 据 ;DEM 与 气象 数据 .CDP 及 人 口 网 格 分 布 数 
据 来 自 中 国 科学 院 资 源 环 境 科学 数据 中 心 ( http:// 
www. resdc. cn/ ) 
2.2 景观 指数 的 选取 与 脆弱 度 指数 的 构建 

景观 格局 脆弱 度 包含 敏感 性 和 适应 性 两 方面 。 
敏感 性 即 在 不 同时 间 和 空间 尺度 上 系统 对 于 干扰 的 
反应 ;适应 性 即 景观 系统 对 外 界 的 干扰 产生 敏感 性 
响应 之 后 , 其 最 终 受 到 干扰 后 的 适应 性 调整 能 
HU., 结合 前 人 研究 ,本 文通 过 景观 敏感 度 指数 、 
景观 适应 度 指数 构建 景观 脆弱 度 评价 模型 。 
2.2.1 景观 敏感 度 指 数 ”景观 敏感 度 指数 是 指 景 
观 在 受 自 然 环 境 变化 或 人 为 活动 干扰 时 的 反应 程 
度 , 表 征 发 生生 态 环 境 问 题 的 难 易 程度 ,其 变化 速率 
越 快 ,敏感 性 越 强 "”。 景 观 敏感 度 指 数 (LS1) 计算 
公式 为 : 


LSI = Y U; x V, (1) 
i=l 


式 中 :nm 为 景观 类 型 数目 ;i 为 景观 类 型 ;U, 为 景观 干 
Jt BE TERR V, 为 景观 易 损 度 指数 。 
(1) 景观 干扰 度 指 数 (U;) 
本 文采 用 具有 代表 性 并 且 与 干扰 密切 相关 的 景 
观 类 型 破碎 度 .分 维 数 倒 数 和 优势 度 3 种 指数 为 基 
础 指标 构建 景观 类 型 干扰 度 指数 ,公式 如 下 : 
U, 2 aFN, + bFD, + cDO, (2) 


IP: EN, 为 破碎 度 指数 ,数值 越 小 表明 破碎 度 越 
小 ;FD; 为 分 维 数 倒数 ,其 值 愈 趋 近 于 1 ,几何 形状 越 
简单 ,表明 受 干 扰 越 多 ,反之 愈 趋 近 于 2 ,结果 越 复 
28, H ER TERR ; DO, 为 优势 度 指数 ,表明 受 景观 类 
型 控制 程度 。 因 此 ,这 3 个 指标 构建 出 来 的 干扰 度 
此 数 能 够 较 好 地 反映 景观 受 干扰 程度 。 借 鉴 前 人 相 
关 的 研究 成 果 , 将 权重 分 别 设 定 为 0. 5,0. 3 和 
0,279. 

(2) 景观 易 损 度 指数 (六 ) 

景观 易 损 度 指数 表示 景观 类 型 与 土地 利用 易 损 
程度 的 有 效 结合 ,借鉴 他 人 研究 成 果 情 ) ,结合 专家 
经 验 和 流域 实际 情况 ,将 景观 类 型 赋 以 相对 权重 , 反 
应 各 景观 类 型 的 易 损 度 :未 利用 地 -7, 林 地 、 草 
地 -5 ,耕地 -3, 建 设 用 地 水 域 -1。 
2.2.2 景观 适应 度 指 数 ”景观 生态 脆弱 度 不 仪 与 
敏感 性 有 关 , 而 且 与 本 身 的 适应 性 密切 相关  。 系 
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统 结构 越 复 杂 且 分 布 均匀 ` 多 样 性 和 功能 越 强 , 则 景 
观 生态 系统 越 稳定 ,对 外 界 干扰 的 适应 能 力 越 强 , 反 
之 亦 然 。 研 究 发 现 斑 块 丰富 度 指数 (RPR) 和 物种 丰 
富 度 存在 很 好 的 正 相 关 '“” ; Simpson 多 样 性 指数 
CSIDI) 在 生态 学 中 应 用 广泛 ,其 与 物种 多 样 性 有 着 
紧密 联系 ;Shannon 均匀 指数 (SHE7) 表示 不 同 景观 
类 型 分 配 的 均匀 程度 ,其 越 高 表明 分 配 越 均匀 ,系统 
越 稳定 。 因 此 ,选取 这 3 个 能 够 表达 生态 意义 且 符 
合 研究 区 情况 的 指数 构建 适应 度 指数 (ZL47) : 
LAI - RPR x SIDI x SHEI 


(3) 
2.2.3 ”景观 脆弱 度 指数 ”结合 景观 敏感 度 指数 
(LSD) 和 景观 适应 度 指数 (ZL47) ,构建 景观 脆弱 度 指 
UVD O>? ,其 值 的 大 小 是 流域 景观 格局 生态 


脆弱 情况 的 定量 表达 : 

LVI = LSI(1 -LA7) (4) 
2.3 驱动 力 指标 定量 
2.3.1 土地 利用 动态 变化 模型 ”本 文采 用 庄 大 方 


等 ”提出 的 土地 利用 程度 综合 指数 (上 ) 作为 人 类 
活动 对 景观 格局 干扰 的 明显 指标 ,具体 公式 如 下 : 


(5) 


L =100x Y (A; x Ci 
i=l 


式 中 :为 栅 格 单元 的 土地 利用 综合 程度 指数 ;4, 为 
栅 格 单元 第 i 级 地 类 分 级 指数 ;C, 为 第 i 级 地 类 占 
网 格 全 部 地 类 面积 的 比重 。 根 据 相关 文献 2 ,将 土 
地 利用 分 为 4 个 级 别 ,分 别 是 未 利用 地 -1, 林 地 、 草 
地 水域 -2, 耕 地 -3 ,建设 用 地 -4。 

2.3.2 气候 要 素 变化 趋势 分 析 ”本 文采 用 Sen 趋 
势 度 法 对 流域 1995—2005 年 .2005 一 2015 年 两 个 
时 间 段 的 温度 及 降水 进行 变化 趋势 计算 ,是 计算 序 
列 的 中 值 , 它 可 以 很 好 地 减少 噪声 的 干扰 ,其 公式 如 
下 : 


X; —X 


B = Median J. visi (6) 


J =} 
式 中 :Median 为 中 位 数 函 数 ;x; 和 x; JS TET 
势 的 序列 数据 ,表示 不 同年 份 的 气温 或 者 降水 ;j IT i 
代表 对 应 的 时 间 即 年 份 ;使 用 6 来 判断 温度 和 降水 
量 在 某 一 时 间 段 内 的 升降 水 平 , 当 其 大 于 0 BEL SR E 
升 趋势 ,反之 呈 下 降 趋 势 。 

2.3.3 其 他 因子 选取 与 定量 处 理 ”在 考虑 土地 利 
用 程度 和 气候 因素 同时 ,本 文 还 选取 了 3 期 对 应 时 
间 段 的 GDP 和 人 口 分 布 数据 进行 重 采样 ,与 其 他 数 


据 分 辨 率 保持 一 致 。 通 过 数字 高 程 模 型 (DEM ) 提 
取 流 域 坡度 \ 坡 向 ,通过 ArcGIS 计算 每 个 机 格 点 距 
城市 .道路 ,河流 的 距离 ,生成 距离 分 布 图 ,限于 篇 幅 
以 2005 年 部 分 数据 为 例 ,如 图 2。 
2.4 随机 森林 分 析 

随机 森林 是 一 种 基于 决策 树 的 机 器 学 习 算法 ， 
每 一 棵 树 都 是 基于 Bootstrap 样本 构建 的 ,然后 对 多 
棵 树 的 预测 结果 进行 组 合 。 每 个 分 类 树 使 用 2/3 的 
数据 来 训练 随机 森林 ,其 余 1/3 的 OOB (out of bag) 
样本 被 保留 ,用 以 进行 模型 验证 。 随 机 森林 通过 对 
指定 变量 组 进行 随机 而 其 他 变量 保持 不 变 , 比 较 由 
此 导致 的 00B error 变化 来 对 变量 的 重要 性 进行 生 
成 。 
2.5 数据 处 理 

景观 脆弱 指数 和 土地 利用 程度 分 布 图 都 是 用 
Fragstats 和 ArcGIS10. 2 选取 99 移动 窗口 计算 得 出 ， 
解决 分 析 尺 度 问 题 。 为 了 方便 气候 变化 与 景观 脆弱 
度 变 化 进行 对 比分 析 , 通过 python. 获得 研究 区 
1995 一 2005 年 与 2005 一 2015 年 两 个 时 间 段 的 气候 
变化 趋势 栅 格 分 布 图 。 对 上 述 所 选取 驱动 力 因子 栅 
格 图 通过 ArcGIS 导出 为 ASCH 格式 ,通过 R 进行 相 
关系 数 计 算 并 完成 随机 森林 建 模 及 精度 检验 ,获得 
驱动 力 的 相对 重要 性 。 


3 ”结果 分 析 


3.1 景观 类 型 变化 特征 

利用 RS 和 GIS 技术 ,对 流域 1995 年 .2005 年 
和 2015 年 6 种 景观 类 型 进行 了 解 译 ,得 到 3 期 景观 
分 布 ( 图 3) ,通过 释 置 分 析 得 到 两 个 阶段 的 景观 类 
型 转移 和 矩阵。 从 表 1 可 以 看 出 ,1995 一 2005 年 流域 
整体 景观 分 布 比 较 稳 定 ,主要 以 耕地 和 林地 为 主 , 但 
两 者 随 着 时 间 推 移 都 出 现 大 幅 减少 ,城乡 用 地 明显 
增加 ,水 域 面 积 有 所 减少 。 其 中 ,耕地 主要 转化 为 城 
乡 用 地 和 林地 ,说 明 这 10 a 内 随 着 经 济 的 发 展 , 城 
镇 化 比较 明显 ,同时 部 分 地 区 开展 了 退耕 还 林政 策 ， 
但 是 整体 来 看 ,林地 人 为 破坏 加 剧 , 主 要 表现 在 林地 
转化 为 耕地 和 草地 ,说 明 毁 林 造 田 , 乱 砍 乱 伐 现象 的 
存在 。 水 域 面积 减少 主要 转化 为 耕地 、 林 地 和 城乡 
用 地 ,说 明 人 为 干扰 的 情况 下 出 现 了 水 土 流失 的 问 
题 。 

从 表 2 可 以 看 出 ,2005 一 2015 年 的 变化 趋势 和 
前 10a 大 体 相 似 , 城镇 化 现象 更 为 明显 , 水 域 和 林 
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Fig.2 Distribution of the driving factors affecting landscape pattern vulnerability in 2005 
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Fig.3 Distribution of the landscape types in the Songhua River Basin in Jilin during the period of 1995 —2015 
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5X1 1995—2005 年 景观 类 型 转移 矩阵 


Tab.1 The transfer matrix of landscape types from 1995 to 2005 /km? 
类 型 耕地 林地 草地 水 域 城乡 用 地 未 利用 土地 合计 
耕地 43 416. 92 41.45 7.80 13.05 141.00 1.19 43 621.41 
林地 122.78 42 801.20 48.72 3.43 25.79 1.65 43 003.57 
草地 4.58 2.18 1225.22 0.91 6.49 4.92 1 244.30 
水 域 23.96 10.82 1.22 1 989.50 5.89 0.26 2 031.65 
城乡 用 地 0.57 1.15 0.12 0.03 4 338.24 0 4 340. 11 
未 利用 土地 8.70 1.83 17.03 0 0.50 727.71 755. TI 
合计 43 577.51 42 858.63 1300.11 2 006. 92 4 517.91 735.73 94 996. 81 
R2 2005—2015 年 景观 类 型 转移 矩阵 
Tab.2 The transfer matrix of landscape types from 2005 to 2015 /km 
类 型 耕地 林地 草地 水 域 城乡 用 地 未 利用 土地 合计 
耕地 43 271.81 24.72 1.33 8.35 269.05 2.25 43 577.51 
林地 84. 69 42 652.43 75.48 3.23 34.99 7.81 42 858.63 
草地 6. 10 2.83 1 280.19 2.85 8.09 0.05 1 300.11 
水 域 29. 60 2.13 2.84 1 955.97 3.82 12.56 2 006.92 
城乡 用 地 0.25 0.48 0.03 0 4 517.15 0 4 517.91 
未 利用 土地 4.44 0.16 0.11 0.03 0.05 730. 94 735.33 
合计 43 396.89 42 682.75 1 359.98 1 970. 43 4 833.15 753.61 94 996.81 


地 转化 为 未 利用 地 明显 增加 ,说明 这 一 时 间 段 内 流域 
受气 候 变 化 影响 ,水土 流失 和 荒漠 化 现象 有 所 体现 。 
3.2 ”流域 气候 要 素 变化 

利用 python 和 GIS 技术 ,计算 1995 一 2005 年 及 
2005—2015 年 两 个 时 间 段 的 气候 要 素 变 化 趋势 (图 


数 绝对 值 均 不 足 0.4, 相 关系 数 较 低 , 且 变 动 相对 较 
大 ,其 原因 是 气候 对 景观 变化 为 长 期 影响 ,因此 更 应 
结合 气候 变化 趋势 图 来 看 ,东部 林地 分 散 的 高 脆弱 
区 发 生 在 温度 升 高 和 降水 减少 的 地 区 ,西北 低 脆 弱 
区 分 布 与 温度 降低 、 降 水 较 高 地 区 分 布 一 致 ,说 明 干 


4)。 从 图 4 可 以 明显 看 出 ,不 同 地 区 气候 要 素 变 化 
具有 一 定 差异 , 且 两 个 时 间 段 内 同一 气候 要 素 变 化 
同样 具有 较 大 差异 性 。 前 10 a, 流 域 东 南部 相对 来 
说 增 温 明 显 , 后 十 年 中 部 增 温 比较 明显 。 前 10 a 
内 ,流域 西南 和 西北 两 个 地 区 雨水 持续 增加 ,东南 部 
有 干旱 趋势 ,后 10 a 内 ,流域 西南 地 区 雨水 依旧 增 
加 ,东南 角 的 降水 开始 增多 ,相对 以 前 ,干旱 趋势 有 
所 改观 ,但 东部 依然 有 部 分 地 区 呈现 出 干旱 趋势 。 
将 对 应 年 份 气温 及 降水 分 布 图 与 景观 脆弱 指数 图 进 
行 相关 分 析 ( 表 3) ,各 期 脆弱 指数 与 气候 值 相 关系 


旱 事件 会 对 景观 脆弱 度 分 布 有 一 定 影响 。 
3.3 ”景观 格局 脆弱 度 及 土地 利用 程度 的 时 空 分 布 
利用 ArcGIS 及 移动 窗口 法 获得 3 个 时 期 脆弱 
度 分 布 ( 图 5) 与 土地 利用 程度 (图 6)。 本 文 根 据 相 
关 文 献 并 结合 松花 江 流 域 自然 地 理 状况 ,将 研究 区 
分 为 5 个 等 级 ,分 别 为 低 脆弱 区 (LV 0. 05) , 较 低 
脆弱 区 (0. 05 < LVI «0. 12), 中 度 脆弱 区 (0. 12 < 
LVI «0. 20) , 较 高 脆弱 区 (0. 20 < LVI «0. 28) ,高 脆 
88 DX (LVI >0. 28) ,并 对 不 同时 期 5 个 分 区 占 总 区 域 
的 百分比 进行 了 统计 。 


表 3 脆弱 度 与 驱动 力 相 关系 数 


Tab.3 Correlation coefficients between ecological vulnerability level and driving factors 


LVI 坡度 坡 向 DEM 利用 度 温度 降水 AH GDP 城市 水 系 道路 
LVI(1995 年 ) | 0.158* 0. 021 * 0. 413 * -0.621* -0.340”  0.389* — 0. 330 -0.237*  0.166*  Á0.014* 0.097* 
LVI(2005 年 ) | 0.161* 0. 022 * 0. 416 * -0.622* .-0.323* .0.342* 一 0. 427 -0.221* 0.178* 0.012 0. 098 
LVI(2015 ^E) | 0.159* 0.021* 0. 418 * -0.627* -0.351” 0.298*  -0.278 -0.282*  0.187* 0.013 0. 102" 


TE LLVI 表示 景观 脆弱 度 指数 ; * 表示 在 P<0.01 水 平 ( 双 侧 ) 上 显著 相关 。 
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Fig.5 Distribution of ecological vulnerability level of the landscapes in the Songhua River Basin in Jilin in 1995, 2005 and 2015 
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6 1995—2015 年 吉林 省 松花 江 流域 土地 利用 度 分 布 
Fig.6 Distribution of land use level in the Songhua River Basin in Jilin in 1995 ,2005 and 2015 


从 空间 分 布 的 角度 来 看 ,松花 江 流域 西北 部 脆 
弱 度 较 低 ,中 部 主要 为 中 度 脆弱 ,高 脆弱 地 区 主要 分 
布 于 西北 边缘 地 区 和 东部 地 区 。 城 乡 用 地 和 耕地 与 
其 他 类 型 相 比较 为 有 序 ,所 以 脆弱 度 较 低 ;中 部 和 东 
部 以 林地 为 主 ,脆弱 度 适中 ;东部 地 区 林地 和 草地 、 
耕地 交错 分 布 ,整体 比较 破碎 ,所 以 呈现 出 了 高 脆弱 
度 ; 西 北 和 中 部 地 区 存在 部 分 难以 利用 的 盐碱地 和 
沼泽 ,呈现 出 局 部 高 脆弱 度 。 结 合 图 5 和 图 6 可 以 
看 出 空间 分 布 的 相关 性 。 同 时 ,结合 图 1 也 体现 出 
脆弱 度 的 分 布 与 流域 地 形 有 一 定 的 一 致 性 。 

从 表 4 可 以 看 出 ,1995 一 2015 年 研究 区 的 脆弱 
度 整体 来 说 相对 比较 稳定 ,但 是 极 低 脆 弱 度 和 极 高 
脆弱 度 的 比例 都 逐年 增加 ,中 等 和 较 低 脆弱 度 比例 
下 降 ,呈现 出 两 极 分 化 的 趋势 。 


表 4 脆弱 分 区 各 年 所 占 比 例 及 其 变化 情况 


Tab.4 Proportions and change of the subareas with 


different ecological vulnerability levels / 96 
a - 较 低 中 度 较 高 AW 
年 份 ag x ^ er Bs di px 
fy (Kiss 脆弱 区 脆弱 区 脆弱 区 脆弱 区 
1995 1.03 34.91 55.03 8.71 0.33 
2005 1.06 34. 85 54.87 8.89 0.33 
2015 1.15 34.78 54.63 9.10 0.34 


A5 随机 森林 效果 评价 


Tab.5 Evaluation of the effects of random forest models 


年 份 MSE RMSE R? 准确 率 
1995 0. 085 0. 292 0. 793 79. 28 
2005 0. 088 0. 297 0. 788 78. 83 
2015 0. 089 0. 298 0. 787 78. 74 


E:MSE 表示 均 方差 ,RMSE 表示 标准 差 ,R 表示 决定 系数 。 


[m^ 


3.4 模型 评价 及 驱动 力 分析 

3 期 景观 脆弱 度 分 布 与 土地 利用 综合 指数 相关 
系数 绝对 值 均 大 于 0.6, 且 与 人 口 和 GDP 也 存在 相 
对 较 强 关系 ,表明 景观 脆弱 度 受 人 为 活动 影响 明显 ， 
同时 气候 因素 及 其 与 河流 远近 对 脆弱 度 的 分 布 有 一 
定 的 影响 ,与 城市 .道路 及 坡度 坡 向 相关 性 较 低 ( 表 
3) 。 单 纯 的 相关 不 能 代表 因果 ,为 了 更 好 地 揭示 驱 
动机 制 ,对 因子 进行 随机 森林 建 模 , 表 5 中 MSE 及 
RMSE 均 较 低 ,R 及 准确 率 相对 较 高 ,模拟 效果 较为 
理想 ,能 够 达到 进一步 研究 的 要 求 。 从 图 7 可 以 看 
出 3 期 驱动 因子 相对 重要 性 的 变化 ,土地 利用 及 
DEM 的 重要 性 与 相关 分 析 中 一 致 , 均 占 有 很 大 比 
例 ,气候 因素 及 GDP 和 人 口 都 占有 一 定 比例 , 且 在 
不 同时 期 起 到 的 作用 略 有 差别 。 河 流 也 有 一 定 的 作 
用 , 且 在 各 个 时 期 相对 稳定 ,与 城市 及 道路 的 距离 ， 
坡度 、 坡 向 这 4 个 因子 在 3 段 时 期 相对 作用 均 较 小 ， 
而 且 驱 动力 随 着 时 间 变 化 整体 作用 相对 较 稳 定 ,但 
GDP 和 人 口 的 作用 相对 气候 因素 有 上 提 趋 势 。 综 
上 说 明 , 以 土地 利用 程度 为 主导 因素 ,结合 人 口 和 
GDP 可 以 很 好 地 作为 人 为 因素 来 分 析 ,将 这 些 驱 动 
力 的 作用 及 分 布 结合 脆弱 度 分 布 状况 来 看 ,流域 景 
观 在 人 为 干预 的 作用 下 , 随 着 土地 利用 程度 的 增加 ， 
建设 用 地 向 四 周 扩展 ,连接 成 片 ,增加 景观 连通 性 ， 
各 类 景观 分 布 也 更 加 均匀 ,结构 更 加 有 序 , 斑 块 边界 
及 形状 也 更 加 规则 ,降低 了 整个 流域 的 破碎 度 ,提高 
了 景观 稳定 性 及 抗 干扰 能 力 ,景观 资源 会 在 土地 利 
用 程度 高 的 地 区 得 到 更 加 有 效 的 使 用 和 整理 ,因此 
这 些 区 域 脆 弱 度 较 低 。 从 脆弱 度 变化 可 以 看 出 , 流 
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图 7 1995—2015 年 流域 景观 脆弱 度 驱 动力 相对 重要 性 


Fig.7 Relative importance of driving factors on landscape vulnerability level in the drainage basin in 1995 ,2005 and 2015 


域 整体 脆弱 度 变 化 不 大 ,下 降 区 域 主要 分 布 于 流域 
土地 利用 程度 较为 发 达 地 区 ,该 区 域 土地 利用 程度 
呈 微 弱 增 加 趋势 ,城镇 化 进程 明显 ,气候 条 件 较为 稳 
定 , 景 观 利 用 愈 发 合理 ;上 升 区 主要 分 布 于 流域 东部 
林 草 地 ,此 区 域 人 口 及 经 济 发 展 水 平均 较 低 , 由 于 人 
为 管理 强度 较 小 ,在 遭受 干旱 迹象 干扰 和 人 为 破坏 
情况 下 ,土地 利用 程度 进一步 降低 ,因此 脆弱 度 上 
升 , 导 致 了 流域 景观 脆弱 度 又 呈现 出 略微 两 极 化 的 
现象 。 


4 结论 与 讨论 


(1) 吉林 省 松花 江 流域 景观 类 型 主要 以 耕地 和 
林地 为 主 ,整体 来 说 景观 分 布 比较 稳定 ,但 随 着 城市 
化 的 推进 ,各 时 间 段 内 不 同 景观 类 型 之 间 存 在 复杂 
的 交汇 ,城市 化 持续 上 升 ,耕地 和 林地 迅速 减少 , 且 
流域 东部 林地 破碎 愈 发 严重 , 中游 水 域 也 呈现 逐年 
减少 的 趋势 ,暴露 出 荒漠 化 及 水 土 流失 问题 。 

(2) 低 脆弱 地 区 主要 位 于 流域 西北 和 中 部 发 达 
城市 ,这 些 地 区 人 为 活动 较 频繁 ,土地 使 用 更 加 有 
FE ,流域 中 部 由 于 升温 和 略微 的 干旱 ,水 土 有 所 流 
失 ,脆弱 度 略 微 增加 ,东部 林地 遭 到 破坏 , DE ER dA 
严重 ,高 脆弱 度 部 分 逐渐 增加 。 西 北部 分 地 区 的 未 
利用 土地 长 期 难以 利用 ,一 直 处 于 高 脆弱 状态 。 整 
个 流域 脆弱 度 分 布 随 着 时 间 变 化 有 两 极 分 化 趋势 ， 
西北 地 区 低 脆 弱 度 稳定 性 持续 升 高 ,东部 林地 脆弱 
度 愈 发 升 高 。 

(3) 松花 江 流域 脆弱 度 分 布 与 变化 受 人 类 活动 
(土地 利用 、 人 口 及 GDP) 影 响 显著 ,其 中 土地 利用 


程度 是 主导 因素 。 整 体 来 讲 , 土 地 利用 程度 越 高 , 斑 
块 破碎 度 越 低 ,连通 性 越 强 ,对 土地 利用 程度 低 的 地 
区 来 说 景观 资源 利用 越 合理 , 抗 干扰 能 力 越 强 , 系 统 
越 稳定 ;脆弱 度 分 布 与 地 形 有 一 定 的 相关 性 ;流域 不 
同时 间 段 气候 变化 趋势 有 所 差异 ,整体 来 看 ,脆弱 度 
变化 与 气候 变化 有 一 定 相关 性 ,如 局 部 地 区 的 干旱 
干扰 ,导致 景观 类 型 改变 及 斑 块 破碎 ,因此 脆弱 度 也 
有 所 升 高 ;脆弱 度 与 距 河流 距离 有 一 定 关 系 ,其 作用 
在 不 同时 间 段 相对 稳定 ;与 城市 及 道路 距离 .坡度 及 
坡 向 等 关系 不 大 。 

(4) 对 于 流域 景观 脆弱 度 的 改善 应 考虑 :适当 
开发 西北 高 脆弱 的 未 利用 地 ,增加 其 利用 ;保护 中 游 
的 水 域 , 加 大 生态 项 目 投入 ,加 修 水 库 ,保护 其 枢纽 
功能 ;保护 东部 的 林 草 地 ,开展 植树 造林 项 目 ,增加 
其 景观 连通 性 ,减少 破碎 化 。 

本 文 以 吉林 省 松花 江 流域 为 研究 区 ,选取 相关 
景观 指数 构建 景观 脆弱 度 ,结合 人 为 地形 及 气象 等 
数据 ,通过 随机 森林 等 方法 分 析 了 景观 格局 和 鹏 弱 
度 时 空 分布 及 驱动 机 制 ,外 界 驱动 力 的 作用 使 景观 
结构 随 着 土地 利用 模式 的 变化 而 变化 ,进而 影响 流 
域 的 生态 服务 价值 和 景观 脆弱 度 , 即 使 在 人 为 干扰 
最 小 的 情况 下 ,景观 功能 也 可 能 因 气 候 和 天 气 条 件 
的 变化 而 发 生 改 变 。 景 观 功能 的 变化 影响 着 生态 系 
统 内 的 能 量 流动 和 养分 循环 ,生态 脆弱 度 也 随 之 发 
生变 化 。 根 据 上 述 分 析 和 实地 调查 结果 ,笔者 认为 ， 
在 过 去 的 20 a 中 ,研究 区 的 景观 格局 脆弱 度 受 到 了 
自然 和 人 为 因素 的 双重 影响 ,其 中 人 为 因素 是 景观 
脆弱 度 较 强 的 驱动 力 , 且 作用 有 上 升 趋势 ,表明 保护 
景观 生态 , 人 类 的 干预 并 不 是 一 个 绝对 消极 的 
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KOO ,科学 合理 的 人 类 保护 措施 可 以 使 景观 格局 更 
加 有 序 , 斑 块 形状 更 加 规则 定 ) ,资源 使 用 更 加 合理 ， 
景观 分 布 更 加 均匀 ,更 有 利于 景观 生态 系统 的 健康 
发 展 。 因 此 ,有 必要 在 自然 继承 的 基础 上 , 辅 以 适当 
的 人 工 干预 ,提高 对 景观 格局 变化 和 脆弱 度 的 适 
应 性 。 


受 研 究 者 经 验 和 知识 水 平等 因素 限制 ,权重 的 


赋值 更 多 的 是 反映 景观 类 型 相对 易 损 程度 ,如 何 采 
用 更 好 的 表达 方法 有 竺 完善。 还 应 进一步 加 强 景观 
指数 的 生态 学 意义 解释 ,驱动 力 对 景观 脆弱 度 的 影 
响 也 具有 阶段 性 和 区 域 性 特征 ,建立 更 完善 的 基于 
生态 过 程 机 理 的 景观 脆弱 度 及 更 完整 时 间 序 列 的 驱 
动力 研究 是 今后 努力 的 方向 。 
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Change of Landscape Pattern Vulnerability in the Songhua River Basin 
in Jilin Province and Its Driving Forces 


SUN Hong-chao, ZHANG Zheng-xiang 
(School of Geographical Science , Northeast Normal University , Changchun 130024 , Jilin , China) 


Abstract: The Songhua River Basin in Jilin Province was taken as the study area , and the three Landsat TM/ETM 
images in 1995 ,2005 and 2015 were used as the data sources. The landscape pattern , spatiotemporal distribution of 
ecological vulnerability and its driving forces were studied. The main conclusions were as follows :(1) Climate change 
in the drainage basin was different from different periods. The dominant landscapes were cultivated land and wood- 
land ,the areas of urban and rural land and grassland were steadily enlarged , but the areas of woodland , cultivated 
land and waters were in a reduction trend. (2) The areas with low and moderate ecological vulnerability were domi- 
nant in the middle and lower reaches of the drainage basin, and the areas with high ecological vulnerability were 
mainly distributed in the southeast of the basin , northwest desertified land in the lower reaches , woodland and grass- 
land. 3) Ecological vulnerability of the drainage basin was relatively stable and in slight increase and polarization 
trends. The distribution of ecological vulnerability was closely related to the human activities, and it was related to 
the topographical conditions to some extent. The effect of climate change was not significant , which , however , could 
still reflect some ecological problems. 

Key words: landscape pattern; vulnerability; spatiotemporal distribution; driving factor; land use; random for- 
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